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Abstract 

The translocation of a macromolecule through a nanometer-sized pore is an interesting pro- 
cess with important applications in the development of biosensors for single-molecule analysis 
and in drug delivery and gene therapy. We have carried out a molecular dynamics simula- 
tion study of electrophoretic translocation of a charged polymer through an artificial nanopore 
to explore the feasibility of semiconductor-based nanopore devices for ultra-fast DNA se- 
quencing. The polymer is represented by a simple bead-spring model designed to yield an 
appropriate coarse-grained description of the phosphate backbone of DNA in salt-free aque- 
ous solution. A detailed analysis of single translocation event is presented to assess whether 
the passage of individual ions through the pore can be detected by a nanoscale field-effect 
transistor by measuring variations in electrostatic potential during polymer translocation. We 
find that it is possible to identify single events corresponding to the passage of counterions 
through the pore, but that discrimination of individual ions on the polymer chain based on 
variations in electrostatic potential is problematic. Several distinct stages in the translocation 
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process are identified, characterized by changes in polymer conformation and by variations in 
the magnitude and direction of the internal electric field induced by the fluctuating charge dis- 
tribution. The dependence of the condensed fraction of counterions on Bjerrum length leads to 
significant changes in polymer conformation, which profoundly affect the dynamics of elec- 
trophoresis and translocation. 
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1 Introduction 



In 1993, Branton and Deamer demonstrated that DNA could be threaded through a biological pore 
and that, by measuring the variations in the accompanying i onic current, information about DNA 
secondary structure could be obtained ( Kasianowicz et all Il996r) . Although these results were 
very promising, little progress has been made towards actual sequencing of DNA due to the limited 
voltage range that can be applied across a biological pore and the difficult y in measuring the cu rrent 



variations, because the shot noise is comparable to the expected s ignal 



alternative approach currently under investigation (Bok hari et al 



(Bokhari et al 



2002). An 



2002) consists in using the new 



capabilities in nanoscale semiconductor technology to create a more optimized pore geometry. A 
device consisting of a nanometer-sized pore in a thin silicon membrane with vertical transistors 
positioned along the wall of the pore may, in principle, be able to detect charge passing through 
the pore by measuring image charges in the transistor. By reducing the noise and maximizing the 
charge and current sensitivity, the proposed geometry may achieve DNA sequencing with order of 
magnitude improvements in speed and cost and minimal pre- and post-measurement procedures 
compared to current sequencing techniques. The expected scientific and therapeutic benefits will 
be dramatic, for example enabling rapid routine screening for treatable genetic disorders. 

Although an empirical approach based on trial and error can lead to some progress in the fabri- 
cation of such devices, a fundamental physical understanding of the delicate interplay between the 
physico-chemical behavior of the macromolecule in solution and the interaction with the confining 
geometry induced by the solid surface of the pore is highly desirable. Molecular modeling is a 
powerful tool to investigate the collective behavior of systems resulting from the complex interac- 
tions between their individual components. Such an approach may allow us to gain insight into the 
phenomena of driven polymer translocation through a nanometer-scale pore from both an applied 
and a fundamental point of view. From an applied point of view, we can assess the feasibility of 
sensing and identifying individual nucleotides, to aid in the design of semiconductor-based de- 
vices and in the interpretation of the experimental measurements. In order to achieve this goal, 
atomistic modeling using high accuracy molecular models must be used to represent the properties 
of single stranded DNA, aqueous solution and substrate quantitatively. With this pre-requisite, it 
is possible to study the local structure, conformation, and short time dynamics of single-stranded 
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DNA in confined geometries, and to investigate the translocation of ions and counterions across 
the channel as a function of the pore diameter and shape and of the magnitude of the potential dif- 
ference applied. From a more fundamental point of view, it is of interest to study the mechanism 
of polymer translocation and the behavior of polyelectrolytes in confined geometries. This is of 
relevance to the process of pore entry, including studies of the role of counterion charge density 
gradients and of the role played by hydrodynamic effects. A better understanding of this process 
can guide the design of pore geometries that effectively steer DNA molecules toward the pore. 

Polymer translocation involves traversal of a free energy barrier arising from several distinct 
physical effects, i ncluding loss of conformat ional entropy of the macromolecule upon traversal (an 



entropic barrier) (Muthukumar 



1999 



2001), mismat ch in dielectric co nstant between the aqueous 



solution and the silicon pore (an electrostatic barrier) (IParsegianL 



1969) an d specific interactions of 



1992J). In this study, 



the macromolecule with the pore surface (an enthalpic barrier) (Zi mm et al. . 
we use a coarse-grained approach, in which groups of atoms are represented by single interaction 
sites, to capture many of the essential physical features of the system. By dramatically decreasing 
the number of degrees of freedom present, this model can engender improved basic understanding 
of field-driven polymer translocation over long length and time scales. 

The remainder of this article is organized as follows: In section II, we define the interaction 
potential used in this study, the parametrization used to represent the phosphate backbone of DNA 
in salt-free aqueous solution, and the simulation methodology. In section III we present results 
obtained on the translocation process and on the influence of Bjerrum length and external electric 
field on the conformational characteristics of the polymer. Finally, conclusions drawn from this 
study as well as future research directions are discussed in section IV. 



2 Methodology 

2.1 Interaction potential 

We study the translocation of a single negatively charged polymer chain designed to represent the 
phosphate backbone of a DNA fragment in an aqueous solution containing sodium counterions. 
Each molecule in this three-species mixture is represented as a set of spherical interaction sites. 
A coarse-grained approach in which groups of atoms are represented by single interaction sites is 
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used in order to reduce the computational cost and to enable the study of the translocation process 
on long timescales. 

We choose a simple interaction potential of the form: U(r N ) = U st r+U v dw+U cou i+U ext +U wai ii. 
The first term decribes intramolecular valence interactions, with U stT representing bond stretching 
interactions. We chose to model the phosphate backbone as a flexible chain without bond angle 
bending or dihedral torsion interactions in our model. It must be pointed out that the chain is not 
fully flexible due to intramolecular van der Waals and Coulombic interactions between second and 
more distant neighbors. The last four terms describe nonbonded interactions, with U v( i w , U cou \, 
U cxt and U wa n representing van der Waals, Coulombic, external electric field, and molecular site 
- substrate interactions respectively. The pore is defined as a cylindrical channel (with axis along 
the z direction) in a solid substrate, which is treated as a smooth van der Waals surface, and the 
external electric field E driving polymer translocation is in the +z direction. These contributions 
to the total potential are defined as: 
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The primes on the sums in Eqs. © and © indicate that nearest-neighbor intramolecular 
nonbonded interactions are excluded from the sums. The distance between sites i and j is 
defined as = | r ^ | = \tj — Tj| where r-j is the position of site i and Z{ its vertical component. 

The remaining step in building the molecular model is fixing the interaction parameters. These 
parameters need to be provided for every distinct combination of site types. In the next subsection, 
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the site types and the interaction parameters used in the case of the coarse-grained model will be 
specified. 



2.2 Parametrization of the coarse-grained model 

An appropriate parametrization of the coarse-grained molecular model is obtained via an approxi- 
mate mapping from the corresponding physical system. The DNA backbone is negatively charged 
and the concentration of positively-charged sodium ions is such that the system is neutral overall. 

Water molecules are treated as point van der Waals (Lennard- Jones) interaction sites, with the 
Lennard- Jones (LJ) parameters determin ed from a 'correspo nding states' mapping of the triple 



point of water onto that of the LJ system djohnson et al. . 



19921) yielding a w = 0.29 nm and e v 



0.79 kcal/mol. These values are taken as the unit of length cr an d the unit of energy €q. The mass 
of a single water molecule is the unit of mass, m = 18 a.u. Using these values, a unit of time 



*o = A/m (To 2 /eo ~ 7 x 10~ 10 s, a unit of charge q = ^eo^o — e/12, and a unit of electric field 
Eq = eo/(qo&o) — 1400 Y/fim are defined, where e is the magnitude of the electron charge. 

The sodium ions are treated as point interaction sites with the same LJ parameters as water 
and with a reduced charge q Na + = aq . We have investigated systems with a = 1.2, 2.4 and 6 
corresponding to a charge q^ a + = e/10, e/5 and e/2. These values have been chosen in order to 
take into account dielectric screening effects in an approximate way. For simplicity, the mass of 
sodium ions is set equal to the mass of water molecules. 

The phosphate backbone consists of alternating phosphate and deoxyribose groups. Each group 
is represented by a van der Waals point interaction site (Figure [T]) with van der Waals parameters 
cr ph = 2cr and e ph = e . The phosphate sites carry reduced electric charge q ph = —q^ a + while 
the deoxyribose sites do not carry any electric charge. The bond stretching parameters are set 
to the following values: r eq = do and k r = 24 kcal/mol. The phosphate chain is quite flexible 
although not fully flexible due to the van der Waals and Coulombic interactions between second 
and more distant neighbors. The masses of both phosphate and deoxyribose sites are set to 5 m . 
Geometrical combination rules for a and e were applied to parametrize the LJ interactions between 
unlike species. 

The pore is defined as a cylindrical channel in a solid substrate, which is treated as a smooth 
van der Waals surface, the atomic sites interacting with the nearest point on the surface of the pore 
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or substrate. For simplicity, the LJ parameters of the wall sites are the same as those of water. 



The diameter of the pore is d 



pore 



5 ctq with a length I 



pore 



2d pore . The diameter of the pore 



is 2.5 times larger than that of the polymer in order to accomodate both the polymer chain and at 
most one solvation sphere. In Branto n's experiments, the alpha-hemolysin pore d iameter has been 



estimated to be between 18 to 26 A dKasianowicz et al 



1996; 



Mel leretal. . 



200 1 ) which is 2 to 2.5 



times the diameter of a single- stranded DNA fragment (D ds = 20A, D ss = 10A dTinland et al 
1997J)). The length of the pore was estimated to be / por e = 52A, in good agreement with the value 
chosen in our simulation. Reduced variables defined in this subsection are used in the remainder 
of this article. 



2.3 Simulation methodology 

We carry out NVT (constant particle number, volume, and temperature) molecular dynamics sim- 
ulations of a single charged polymer in solution, under an external field, with periodic boundary 
conditions. The unit cell contains a reservoir separated by an impenetrable substrate of thickness 
/p 0re and surface area equal to the lateral dimension of the unit cell, and containing a cylindrical 
hole (nanopore) of diameter G? pore (Figure EJ- 

The chain consists of 40 sites, alternating between phosphate and deoxyribose groups, and the 
initial state is built such that the charged polymer is in the reservoir outside the pore cavity. The 
system contains, in addition to the polymer, 6000 water molecules and 20 Na + counterions with 
initial positions chosen randomly in the unit cell. The simulation was carried out at a reduced 
density p* = 0.85 and a reduced temperature T* = ksT '/e = 0.75, corresponding to the dense 
fluid phase of a LJ system. The applied external electric field of magnitude E* = E/E = 0.5, 
corresponding to a field E =709 V//im, is turned on after an initial equilibration run of duration 
t* = 500. In the experiments conducted with the alpha - hemolysin nanopore, a voltage of 70-300 
mV was applied across the system. The electrodes are located centimeters away from the pore, but 
the bulk of the voltage drop occurs across the biological membrane, leading to an estimated electric 
field of 15-60 V/fim, at least an order of magnitude smaller than the field used in the present work. 
Such a large field is used in the present study in order to ensure that field-driven translocation 
occurs on a time scale accessible to simulations. 

A single-timestep velocity- Verlet MD integration scheme was used to integrate the equations 
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of m otion with a timestep St = 2.5 x 10 3 , and the weak-coupling algorithm of Berendsen et 



al. (Berendsen et al 



1984) was used to maintain constant temperature. All LJ interactions were 



truncated at 2.5 a and no long-range corrections were applied. Long-range Coulomb interactions 



are evaluated to high 



Essman 


n et a . 


jEwald 


1921: 



accuracy using the particle-mesh Ewald (PME) method (Dar den et al 



1993; 



1995 ). The PME tech nique, which is derived from the conventional Ewald method 



Toukmaji et al. 



1996), makes use of the EFT to efficiently compute the long range 



part of the electrostatic interaction. In this study we used a highly optimized version of the PME 
with a relative accuracy of 10~ 4 . 

We have carried out 3 simulations of systems with reduced charge magnitude q* = 1.2, 2.4 and 
6.0 for durations of 4.2 x 10 6 , 3.0 x 10 6 and 5.2 x 10 6 timesteps, respectively. 



3 Results 

The translocation process consists of three distinct stages: In the first stage, the polymer drifts 
through the reservoir, driven by the external electric field. In the second stage, the polymer is 
pushed against the lower substrate surface and behaves like a 2-dimensional chain wandering in 
search of the pore entrance. In the third stage, the polymer finds the entrance to the pore and 
translocation takes place. 

Due to the relatively small system size, the time spent by the polymer in the second stage is 
small compared to that in the two other stages. We first study the electrophoresis of the charged 
polymer in the reservoir, then the kinetics of the translocation process and the possibility of iden- 
tifying individual ions passing through the pore from the electrostatic signal produced. 



3.1 Electrophoresis in the reservoir 



Electr ophoresis is a complex process dependent on the electrophoretic friction coefficient ( Ma nning , 
1981). This coefficient arises both from the macroion/solvent interactions, comprising the monomer/solvent 
interaction (the intrinsic friction) and from hydrodynamic interactions between monomer pairs, 
and from macroion/ion interactions. Sinc e the macroion is ne gatively charged, it is surrounded by 



an oppositely charged liquid atmosphere (Fu oss et al. . 



1959). The force acting on the ions in the 



atmosphere is transmitted to the solvent, and therefore the atmospheric liquid behaves as a charged 
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volume. Under an applied electric field, this charged liquid is subject to a volume force in the 
opposite direction from the drift velocity of the macroion. The macroion thus moves against a 
local hydrodynamic flow that slows its motion (relative to a hypothetical drift without the charged 
liquid). In addition, the electric field deforms the charged volume, increasing the atmospheric 
charge density at the end of the polymer opposite to the end which drift and decreasing ahead of 
the macroion. This effect induces an internal field (often referred in the literature as the relaxation 
field, or more accurately, the asymmetry field) opposing the external electric field, which also 
slows down the macroion motion. 

Counterion condensation is in large part responsible for the structure and behavior of the liquid 
atmosphere and , as r ecognized in the late 60's independently by Oosawa and Manning ( Oosawa . 



197 It Manning , 



19781) . plays a crucial role in the strong attractive interaction that acts between 



highly charged macroions such as DNA, leading to significant conformational changes. Counterion 
condensation results in a competition between Coulombic energy and entropy in minimizing the 
free energy of an aqueous solution containing mobile ions in the vicinity of an isolated macroion. 
For rodlike objects, whether the entropy or the Coulombic attraction dominates depends on the 
charge density of the macroion measured in terms of the Bjerrum length \b = q 2 /ksT, defined as 
the distance at which the Coulombic interaction between two charges is equal to the thermal energy 
ksT. Counterion condensation occurs when the distance b between charges on the macroion is 
small enough for the dimensionless parameter Y = \ B /b (the Manning parameter) to exceed unity. 
With the reduced units introduced in the previous section, a reduced Bjerrum length A^ = q* 2 /T* 
is defined. In the present model, the distance between 2 charges is b = 2a , leading to a Manning 
parameter T = q* 2 /2T*. Values for A^ and T are listed in Table [T]for the systems with reduced 
charge q* = 1.2, 2.4, 6.0. At the reduced temperature chosen in this study, Coulombic interactions 
dominate over the thermal energy for q* = 6.0 and q* = 2.4 systems, which thus should exhibit 
strong counterion condensation. For the q* = 1.2 system, T = 0.96, which is close enough to 
unity to also exhibit counterion condensation. However, in this case we expect to have a much less 
strongly bound counterion - macroion complex than for q* = 6.0 and q* = 2.4. 

It is difficult to d efine an unambiguous criterion for condensation. Indeed, as pointed by Oo- 



sawa (Oosawa, 



19711) . for a coiled polymer chain each charged group makes a sharp and deep 



potential hole at its position, each linear part of the chain makes a sharp and deep potential valley 
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along its length and the coiled chain as a whole makes a potential trough in its apparent volume. 
Counterions located in these 3 regions are considered to be condensed, with the counterions at 
charged group holes being localized and the counterions in the other two regions being mobile. 
Therefore, an estimate of the degree of condensation depends on the criterion used to discriminate 
between free and condensed counterions. In the present study, a counterion is said to be condensed 
if it is within a distance r* from an ion of the polymer. 

Figure |3] shows representative configurations for the three systems both with and without ap- 
plied external electric field E*. We first study the system behavior without an external applied 
field. The three systems exhibit counterion condensation but the nature of the resulting macroion- 
counterion complex varies with the magnitude of the charge. For q* = 6.0, a large majority of the 
counterions are closely associated with the polymer chain, being mostly bound or in the narrow 
potential valley along the length of the chain, while for q* = 1.2 the counterions are mostly lo- 
cated in the potential valley of the polymer chain. An intermediate situation occurs for q* = 2.4. 
Figure |4] shows the radial pair correlation function gi C (r*) between ions on the chain and counte- 
rions, for q* = 1.2, 2.4 and 6.0 and E* = 0. The common features exhibited are a main peak 
at r* ~ 1.5, corresponding to localized condensed counterions, and a shoulder around r* ~ 2.0, 
corresponding to mobile counterions in the narrow potential valley around the chain. The height of 
the main peak for q* = 6.0 is roughly 7 times larger than for q* = 1.2 and 4 times larger than for 
q* = 2.4, confirming a stronger counterion localization for larger charge magnitude. For q* = 1.2, 
Figure |5] shows the relative importance of longer range correlations with respect to short range 
correlations with the presence of secondary peaks around r ~ 2.3<7o and r ~ 3.5a , indicating that 
a significant counterion fraction is located in the potential trough of the chain. Table[T]displays the 
average number of condensed counterions for the 3 systems for different values of r* correspond- 
ing roughly to the bottom of the main peak in g ic (r* = 1.8), the shoulder (r* = 2.2) and a value 
taking into account the long range correlations (r* = 6.0). 

Under the applied external field, the counterions remain located very close to the macroion 
forming bound dipoles for q* = 6.0. In addition, a conformational change occurs with the forma- 
tion of a collapsed polymer structure. By contrast, for q* = 1.2, the counterions are not condensed 
anymore (Figure 0). The force produced by the applied external field overcomes the attractive 
Coulombic force and strips the counterions from the chain. As a first approximation, the reduced 
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critical charge at which the Coulombic force exactly balances the force produced by the exter- 
nal field is q* = E*f* 2 where f* is the average reduced distance between a polymer ion and a 
counterion. Assuming a condensed state with all counterions bound to the polymer, the estimate 
f* — 1.5 - 2.0 leads to q* = 1.125 - 2.0. Thus, no counterions will be stripped from the chain 
for q* = 6.0, while a large counterion fraction will be stripped away for q* = 1.2. In the latter 
case, q* is a lower limit since the counterions are not closely bound to the chain (see Figure 0). A 
Monte Carlo simulation study of the counterion con densation on a spher ical macroion has shown 



the same phenomenon under a strong external field (Tan aka et al. 



2002). With this estimate, the 



counterion fraction located in the potential trough of the chain will be stripped from the polymer 
for q* = 2.4. 

Figure |5] shows the ion-counterion pair correlation function when an external electric field is 
applied on the systems. For q* = 1.2, the magnitude of the main peak for q* = 1.2 is reduced 
by an order of magnitude, confirming that a large counterion fraction is stripped from the chain. 
The electric field has a much less drastic effect on the two other systems: A small fraction of 
counterions are stripped from the chain for q* = 2.4 while for q* = 6.0, the field enhances slightly 
the counterion binding to the polymer. Average values for the number of condensed counterions 
under external electric field are given in Table [T] 

The main qualitative features exhibited during the time evolution are independent of the chosen 
value for r*. Figure |^1 shows the time evolution of the number of condensed counterions for the 
three systems for r* = 1.8 and r* = 6.0. The r* = 6.0 system does not exhibit any significant 
changes in the number of condensed counterions over the duration of the simulation. The coun- 
terions are mostly condensed on the chain, and the electric field slightly enhances this effect. By 
contrast, for the q* = 1.2 system, a strong decrease in the number of condensed counterions occurs 
in an external field, the new equilibrium configuration being reached after a time t* ~ 2500. The 
q* = 2.4 system exhibits an initial decrease of the number of counterions induced by the external 
electric field similar to that for q* = 1.2 with a new equilibrium state also reached after a duration 
t* ~ 2500. This effect is weaker, however (see Tableland Figure |5} due to the fact that a larger 
fraction of counterions are condensed on the chain (q* < q* = 2.4). The behavior during the 
translocation process is discussed in the next subsection. 

We turn now to the study of the polymer conformations as a function of magnitude of the site 
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charge q* and electric field E*. A global measure of the polymer conformation is the mean square 
radius of gyration R 2 , defined as: 



Rz 



1 N 

1=1 



< (n - RY > 



(6) 



where is the position of the i th polymer site, N is the number of sites in the chain, and R is 
the position of the center of mass of the chain. Average values of R 2 g over a duration t* = 2500 
for the 3 systems with and without applied electric field are listed in Table [H Figure |7] shows the 



time evolution of R 2 for q* = 1.2, 2.4 and 6.0. For q* = 1.2, the chain becomes more elongated 
under the action of the external field because the mobile counterions are stripped away by the 
field, leaving the bare negative ions on the polymer and leading to an electrostatic stiffening of the 
polymer. For q* = 2.4, the first effect of the electric field, before translocation through the pore, is 
to strip away the fraction of counterions which are in the potential trough of the polymer, resulting 
in a loss of electrostatic screening of the negative charges on the polymer and an increase in R 2 . 
Then the polymer chain is driven through the reservoir by the electric field and pushed against the 



lower substrate (FigureOJ). In this regime, the p olymer behaves 



1997 



ike a 2-dimensional chain and its 



1998). Conformational changes during 



radius of gyration may be significantly reduced (IFooL 
the translocation process are discussed in the next subsection. 

For q* = 6.0, the radius of gyration is of the same order of magnitude as for q* = 1.2 and 
q* = 2.4 when no external field is applied. Under an external electric field, no significant change 
in the magnitude of the radius of gyration is observed over a time interval t* ~ 5000, then a 
strong decrease in R 2 occurs, indicative of the tightly collapsed structure adopted by the chain 
(Figure 13). Additional studies are needed to determine whether the observed collapse is induced 
by the external field or is simply caused by strong counterion condensation. Biological processes 



and theoretical works tend to favor the 



packing of DNA into a cell (Bloomfield, 



atter h ypothesis. The chain collapse is analogous to the 



1996). This packing requires overcoming an enormous 



Coulombic barrier in a highly dilute aqueous solution containing a small concentration of polvalent 
cations. Simulations and theo ry have demonstrated short range attraction between tw o macroions 



mode 



ed as charged cylinders (iGronbech- Jensen et al 



mers ( St evens , 



1997; 



Lvubartsev et al. 



1995) or stiff poly- 



1999). This attraction has been ascribec 



induced by counterion condensation (H a et al 



1997; 



to correlated fluctuations of the counterions 



Manning et al. 



1994). Self-attraction have 
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been shown to occur also in flexible polyelectrolytes ( Stevens et al 



Brilliantov et. al 



1995; 



Schiessel et al 



1998 



1998). Theoretical work using simple scaling arguments have shown that coun- 



terion condensation modifies the second virial coefficient of a polyelectrolyte due to the fact that 
at low enough temperature (or at high enough counterion valence) the counterions approach close 
enough to the macroio n to form dipoles, leading to charge - dipole and dipole - d ipole interactions 



(Schi essel et al 



1998). Similar ideas have been developed by Brilliantov et al. (IBrilliantov et al 



1998), who also predict also a first order phase transition between a stretched polyelectrolyte and a 
strongly collapsed polyelectrolyte as a function of magnitude of electric charges due to counterion 
condensation. 



3.2 Translocation process 



The thermodynamics of electrophoresis in th e presence of a narrow pore is complex and involves 



traversal of entropic and enthalpic barriers 



Slonkina et a 



Kumar et al., 



2002; 



200o[ 



Lubenskv et al. 



Lee et al. 



1999; 



(Muthuku 



mar 



Sung et al 



1999 



1996; 



2001 



Boehm, 



Ambiornson et al 



199 



Sebastian et al. 



2002; 



2000; 



2001). When the chain enters the pore, its conformational entropy 



is reduced, due to chain elongation. On the other hand, due to the small size of the pore aperture a 
significant fraction of the counterions are stripped from the chain and gain configurational entropy. 
The counterion unbinding also results in an increase in the Coulombic energy. In addition, the elec- 
trostatic energy of interaction of the macroion with the external field decreases as the negatively 
charged chain moves through the pore, and the electrostatic interaction energy of the free coun- 
terions also decreases. The situation is reversed when the polymer exits the pore. However, the 
gain in conformational entropy of the polymer upon exiting the pore is probably greater than the 
entropy loss associated with pore entry, since before pore entry the polymer is pushed against the 
lower substrate surface (z = —Zq) and behaves like a conformationally restricted two-dimensional 
chain while it behaves like a three-dimensional chain upon exiting the pore. Detailed free energy 
calculations are needed to study the importance of these contributions in the translocation process. 
We focus here on the kinetics of the translocation process under a large electric field. 

Translocation was only observed in the q* = 2.4 simulation. As we have seen in the previous 
subsection, for q* = 6.0, strong counterion condensation results in repulsive screened Coulomb 
interactions between the ions on the chain and self-attraction due to charge-dipole, dipole-dipole 
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interactions and/or charge fluctuation along the chain. In the q* = 6.0 system, the chain adopts a 
collapsed structure with an effective diameter larger than the pore diameter, effectively preventing 
translocation. There is, a priori, no reason for the q* = 1.2 system not to translocate since the chain 
adopts an extended conformation. The simulation has to be performed over a long enough period 
of time to permit the polymer to find the pore entrance via diffusion along the substrate surface. 

Due to the large external field, translocation of the polymer in the q* = 2.4 system occurs after 
a short drift time spent in the reservoir. The polymer finds the entrance of the pore after a time 
interval At* = 4000 (At = 2.8 ns) after application of the field. The upper part of Figure[5] shows 
several stages in the translocation process. For clarity, only the polymer chain and the counterions 
are shown. The translocation starts when the first bead (head) enters the pore at z = —z and ends 
when the last bead (tail) exits the pore at z = z . The total translocation process occurs over a 
period of time At* ~ 950, corresponding to 0.665 ns. Timescales reported for translocation of 
ions and water molecules are in the range 10~ 6 — 10~ 9 s, in reasonable agreement with our results. 
However this time is much shorter than the characteristic time reported for ssD NA translocation 



(10 6 s) (Kasian owicz et al 



1996; M eller et al. 



2001; 



Henrickson et al. 



2000) due to the larger 



electric field used in the present work. 



2002J) in 



We have decomposed the translocation process into three stages ( Slonkina et al. 
which the conformation of the polymer and its interactions with the pore are significantly dif- 
ferent: the first stage corresponds to the translocation of the head from the entrance to the exit 
of the pore, the second stage corresponds to the exit of the head from the pore and the translo- 
cation of the tail from the reservoir to the entrance of the pore, and the third stage corresponds 
to the translocation of the tail from the entrance to the exit of the pore (see Figure [5]). The du- 
ration of each stage of translocation is: At* = t^ - t^ zn = 145.74 (At t = 0.102 ns); 



a*2 = 



Ahead 



424.75 (At 2 = 0.297 ns) and Att 



=-20 

j-tail j- 

''2=20 1 



tail 



379 (At 3 = 0.266 



ns). The duration of the second and third stages are comparable, but are nearly three times as long 
as the first stage. These differences in time are probably due to variation in the magnitude of the 
total electric field i^otai = E* + E* nt experienced by the polymer during the translocation, where 
E* nt is the internal field induced by the charge distribution. 

The bottom part of Figure [S] shows the reduced electrostatic potential V* = V/V (V = e / qo) 
due to ions and counterions in the (xz) plane passing through the middle of the pore. We clearly 
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see an inversion in the sign of the potential difference across the pore during the second stage of 
polymer translocation due to a change in the relative charge density on the two substrate surfaces. 
This change leads to an inversion in the direction of the internal field E* nt , which in turn slows 
down the translocation process. The upper and lower substrate surfaces behave like a capacitor, 
with an excess of positive charges on the upper surface - the free counterions trying to translocate 
and pushed against the surface by the external electric field - and an excess of negative charges on 
the lower surface - the polymer chain wandering at the pore entrance. This situation is reversed 
when enough negatively charged monomers have translocated, analogous to the discharge of a 
capacitor. 

Figure |9] shows quantitatively the time evolution of the reduced internal electric field E* nt , 
defined as E* nt = AV^ ore /l* ore , across the pore during the translocation event, with AV^* orc being 
the average potential difference across the pore defined as Ay p * ore = 1/4 

V*(x i} y u z )) with (x { = 0; y { = ±r ) for i = 1,2 and (x { = ±r ; y { = 0) for i = 3,4 (see 
Figure |2j). It confirms that the internal field decreases in magnitude and changes sign when the 
middle of the chain crosses the middle of the pore (i.e., in the second stage) resulting in a smaller 
magnitude of total electric field E* ot and a slower translocation process. 

We turn now to the identification of single translocation events involving either a counterion 
or charged sites on the polymer. In order to mimic the signal recorded at the gate of a transistor 
located on the inner cylindrical surface of the pore, we have computed the electrostatic potential 
V* at 4 equidistant positions (x — 0; y — ±r ) and (x = ±r ; y — 0) on the surface of the 
mid-plane (z = 0) of the pore. We have simultaneously recorded events corresponding to the 
passage of charged sites (either counterions or charged polymer beads) through the z = plane 
to investigate the correlations between signal and site positions. Figures and [TT] show the time 
evolution of V* averaged over the four sensors as well as ion passage events during the transloca- 
tion process. A positive spike records the passage of a charge in the same direction as the external 
electric field while a negative spike records the passage of a charge in the opposite direction. No 
counterions are pulled into the pore with the polymer so we do not have any counterions translocat- 
ing against the electric field direction during polymer translocation. A strong correlation is found 
between counterion positions and the signal (sharp peaks in V*). It is also fairly easy to detect 
when the polymer translocation takes place (and its approximate duration), but the discrimination 
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of individual charges on the polymer chain is more problematic due mainly to the small separation 
between neighboring charged sites on the chain, fluctuations in chain position, and the simultane- 
ous translocation of counterions. Further analysis or filtering of the signal is required to achieve 
single-ion discrimination even with this simple model. 

We notice that a significant larger number of counterions find the entrance of the pore during the 
third stage of polymer translocation. During the translocation process, the number of condensed 
counterions on the chain increases significantly (Figure |6j). Due to the applied electric field, an 
excess counterion density is located near the upper substrate surface and a significant fraction of 
them - the ones near the aperture - condense back onto the polymer as soon as the head of the 
chain exits the nanopore (see Figure [8]). The condensed counterions near the tail of the polymer 
are relatively mobile and have a significant probability of being stripped from the chain due to 
the external field and friction with the water molecules. The translocation of these counterions is 
facilitated by the polymer, which acts as a guide to the pore entrance. 

In addition, the head of the polymer exiting the nanopore during the translocation process 
experiences a stronger electric field than inside the pore (see above). The polymer tail is still 
pushed by the field against the lower substrate surface creating an anchoring point on the substrate 
(see Figure [U) while the head of the chain is stretched by the electric field at the pore exit. This 
effect results in an increase in R 2 g (Figure 13). After the head has exited the pore, counterion 
condensation back on the head, screening the Coulombic interactions and reducing the stretching 
force, as well as a weaker anchoring at the entrance of the pore due to the presence of fewer 
monomers, are responsible of the decrease in polymer extension, the chain behaving like a rubber 
band. As the translocation proceeds, the tail enters the pore and the release of the anchoring point 
induces a significant contraction of the chain. After the translocation, the chain behaves as a three- 
dimensional chain, resulting in an increase in R 2 g . 

4 Conclusions 

Using a coarse-grained bead-spring model, we have studied the electrophoretic translocation of 
a charged polymer through a nanopore for various charge magnitudes under the application of an 
external electric field. Three regimes can be identified, corresponding to drift of the polymer in 
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the reservoir, diffusion of the polymer along the substrate surface in search of the pore entrance 
and polymer translocation through the pore. The three systems studied, with reduced charges, 
q* = 1.2, 2.4 and 6.0, respectively, exhibit counterion condensation when no electric field is ap- 
plied, but the q* = 1.2 system presents a much less strongly bound macroion-counterion complex 
than the others two. Counterion condensation leads to significant changes in polymer confoma- 
tion, with a rather extended chain conformation for q* = 1.2 and a collapsed chain conformation 
for q* = 6.0, due to self-attraction induced by charge-dipole/dipole-dipole interactions and/or 
counterion fluctuations along the chain. An external electric field can strip the counterions from 
the chain, depending on the relative magnitude of the field and the charge of the ions/counterions, 
inducing significant changes in polymer conformation. Due to the strong external electric field 
used, translocation occurs quickly for the q* = 2.4 system and simulation over longer time scales 
should lead to the translocation of the q* = 1.2 system. By contrast, translocation in the q* = 6 
system is inhibited due to the collapsed conformation adopted by the chain. The kinetics of the 
polymer translocation are not constant over the whole process and the chain translocates more 
slowly when its midpoint crosses the mid-plane of the pore due to an inversion in the direction 
of the internal electrostatic field produced by the time-varying distribution of charges. The upper 
and lower substrate surfaces behave qualitatively like a capacitor which undergoes a discharge. By 
recording the electrostatic potential we were able to identify the passage of individual counterions 
unambiguously. However, it is difficult to identify translocation events involving single charged 
sites on the polymer due mainly to the small distance between neighboring ions, fluctuations in 
polymer position, and simultaneous counterion translocation. 

This work is mostly exploratory and many questions remain unanswered. Future work will 
focus on both the translocation phenomenon and on bulk electrophoresis. Using the simple model 
presented in this work, we will study the electrophoresis process in detail, including the effect of 
external electric field on counterion condensation as a function of the ion charge magnitude and 
macroion concentration. Special emphasis will be given to the modification of macroion conforma- 
tion induced by the external field. In the framework of the simple model introduced in this study, 
we will study the translocation process statistically by steering the polymer through the pore and 
performing free energy computations to assess the relative importance of entropic and enthalpic 
contributions to the free energy of pore traversal. Since one of our aims is to study the feasibility of 
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a biosensor for fast DNA sequencing we have initiated large-scale atomistic simulations of ssDNA 
in aqueous solution to study the electrostatic field produced by different nucleotides in the vicinity 
of the macromolecule, to explore the possibility of discriminating different bases based on their 
electrostatic signature. 
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Figure 1: Mapping of the phosphate backbone (left) onto a coarse-grained model. The phosphate 
group is represented by a spherical interaction site (red) carrying a charge q = —q^ Si + and the 
deoxyribose group is represented by a spherical interaction site (green) carrying no charge. 

Figure 2: Side view (left) and top view (right) of the simulated system, comprised of a reservoir 
with a substrate containing a nanopore of radius r and length 2z . 

Figure 3: Representative configurations of the 3 polyelectrolyte systems, without external field 
(top) at t* = 500 and with applied electric field (bottom) at t* = 10500 for q* = 1.2 and q* = 6.0 
and at t* = 4000 (before translocation) for q* = 2.4. For clarity, the water molecules have not 
been displayed. 
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Figure 4: Ion-counterion pair distribution function gi C (r) for q* = 1.2 (solid line), 2.4 (dashed 
line) and 6.0 (dot-dashed line) without external field. 

Figure 5: Ion-counterion pair distribution function <7i c (r) for q* = 1.2, 2 .4 and 6.0 (from top to 
bottom) for = (solid line) and E* = 0.5 (dashed line). 

Figure 6: Time evolution of the number of condensed counterions for q* = 1.2, 2.4 and 6.0 (from 
top to bottom) using a condensation criterion r* = 1.8 (thick solid line) and r* = 6.0 (thin solid 
line). For q* = 2.4, the arrows represent the time location of the four translocation stages (A, B, 
C, D) presented on Figure [SJ 

Figure 7: Time evolution of the mean square radius of gyration < R 2 g > for q* = 1.2, 2 .4 and 
6.0 (top to bottom). For q* = 2.4, the arrows represent the time location of the four translocation 
stages (A, B, C, D) presented on Figure 

Figure 8: Top: Four stages in the translocation process, from left to right: head of the polymer 
at the entrance of the pore (t* = 4472.75), head at the exit of the pore (t* = 4618.5), tail at 
the entrance (t* = 5043.25), tail at the exit (t* = 5422.25). Bottom: Corresponding reduced 
electrostatic potential V* = V/V (V = e /q ) produced by the charge distribution in the plane 
(xz) passing through the center of the pore. Values of the potential outside the range chosen (-5,+5) 
are assigned to the colors corresponding to the minimum and maximum values, and for clarity the 
water molecules have not been represented. 

Figure 9: Time evolution of the induced electric field E* produced by the distribution of charges. 
The arrows represent the time location of the four translocation stages presented on Figure 

Figure 10: Time evolution of the reduced electrostatic potential V* and of the events corresponding 
to the passage of counterions (top) and ions (bottom) through the midplane of the pore (z = 0), 
during polymer translocation. A positive peak corresponds to an ion or counterion passing in the 
same direction as the field while a negative peak corresponds to an ion or counterion passing in 
opposite direction. 

Figure 1 1 : Time evolution of the reduced electrostatic potential V* and of the events corresponding 
to the passage of ions at z = 0, during the polymer translocation. Magnification of Figure fTOl 
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Table 1: Electrostatic and conformational properties of the simulated systems. 





q* = 1.2 


q* = 2.4 


q* = 6.0 




1.92 


7.68 


48.0 


r 


0.96 


3.8 


24.0 


Number of 
condensed counterions 

(E = 0) 


2.347 ± 0.158 (r* = 1.8) 


4.761 ±0.712(1.8) 


11.785 ± 1.897 (1.8) 


2.897 ± 0.263 (r* = 2.0) 


5.706 ± 1.074 (2.0) 


12.726 ± 2.044 (2.0) 


10.174 ± 2.266 (r* = 6.0) 


12.358 ± 2.428 (6.0) 


16.866 ± 1.648 (6.0) 


Number of 
condensed counterions 
(E* = 0.5) 


0.283 ± 0.138 (r* = 1.8) 


3.029 ±0.902 (1.8) 


14.254 ±0.165 (1.8) 


0.268 ±0.149 (r* = 2.0) 


3.623 ± 1.265 (2.0) 


14.510 ±0.079 (2.0) 


1.567 ±0.494 (r* = 6.0) 


10.404 ±3.158 (6.0) 


18.026 ±0.049 (6.0) 


< R] > (E* = 0) 


23.030 ± 1.375 


25.789 ± 1.459 


20.202 ± 1.008 


<R 2 g > (E* = 0.5) 


60.007 ± 2.716 


60.837 ±7.157 


13.961 ±0.279 
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